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^ Abstract 

We discuss the three-body properties of identical bosons exhibiting large scattering length in 
d two spatial dimensions. Within an effective field theory for resonant interactions, we calculate the 

^> leading non-universal corrections from the two-body effective range to bound-state and scattering 

observables. In particular, we compute the three-body binding energies, the boson-dimer scattering 

^— I properties, and the three-body recombination rate for finite energies. We find significant effective 

range effects for three-body observables in the vicinity of the unitary limit. The implications of 
^ this result for future experiments are briefly discussed. 
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I. INTRODUCTION 



In the last few years, ultracold quantum gases have become a versatile tool to investigate 
few- and many-body phenomena in strongly interacting quantum systems. With the help of 
Feshbach resonances, the atomic interaction strength can be changed at will. This allows, 
for example, for the investigation of the BCS-BEC crossover in Fermi gases, the creation 
of molecules, the observation of the Efimov effect and of other universal phenomena (see, 
e.g., [U 12] and references therein). The possibility of using optical lattices makes them also 
interesting for the simulation of condensed matter problems such as the Hubbard model [3] . 
Special trap geometries allow for the creation of lower-dimensional systems. They can, for 
example, help to understand high-temperature superconductivity which is a two-dimensional 
{2d) problem. Moreover, 2d systems are interesting on their own since their behavior can 
be qualitatively different from three dimensions. 

Here, we concentrate on few-body phenomena in ultracold atomic gases with large scat- 
tering length. Such phenomena are of great interest because they are insensitive to the 
details of the interaction at short distances. They can be described in an expansion around 
the unitary limit. This limit refers to an idealized system where the range of the interaction 
is taken to zero and the scattering length a is infinite. To leading order in this expansion, 
the low-energy observables are universal. They are determined by the scattering length a 
of the particles alone. The leading non-universal corrections are due to the effective range 
of the interaction. In the current paper, we focus on these corrections. Since there is no 
Efimov effect [1] in two dimensions E], three-body interactions are suppressed and enter 
only at higher orders. 

The definition of the scattering length a in 2d is ambiguous since cot 6 diverges loga- 
rithmically as the wave number k approaches zero and different conventions are used in 
the literature. We follow the conventions of Verhaar et al. [7j, in which the effective range 
expansion of the scattering phase shift is given by 

cot 6ik) = l[lE + In (y ) } + + ^(^') ' (1) 

where 7^ ~ 0.577216 is Euler's constant. Note that the scattering length in two dimensions 
is always positive. 

In the limit a ^ |r|, the binding energy of the shallow dimer is universal, 

i?2 = 4e-2^"^ + 0(rV«')- (2) 

The binding energies of three- and four-body states in this limit have been calculated by 
various groups and are also universal |21 El I5HTT]. Since there is no other parameter in the 
problem, the energies must be multiples of the dimer energy. There are two three-body 
bound states which were first calculated by Bruch and Tjon [5]. Their binding energies 
are |H] 

= 1.2704091(1)^2 and = 16.522688(1) , (3) 

where the number in parentheses indicates the numerical error in the last quoted digit. The 
first calculation of the four-body bound states in 2d was carried out by Platter et al. [H]. 
They found also two universal bound states with binding energies 

= 25.5(1)^2 and = 197.3(1) . (4) 
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FIG. 1: (Color online) Spectrum of universal two-, three-, and four-body states in two spatial 
dimensions as a function of the inverse scattering length 1 ja. 



These results were later confirmed in Ref. |T^. In Fig. [T| we illustrate the scattering length 
dependence of this spectrum. The universal few-body states do not cross the continuum 
threshold = for any finite value of the scattering length. In contrast to three dimensions, 
the Id universal states can therefore not be observed as zero-energy resonances in few-body 
recombination. 

For large values of ^ 1, one can derive the universal properties of shallow A^-boson 
ground states close to the unitary limit [S]. In particular, the binding energy of the 
A^-boson ground state increases geometrically with A^: 

- 8.567, A^ > 1. (5) 
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Thus, the separation energy for one particle is approximately 88% of the total binding en- 
ergy. This is in contrast to most other physical systems, where the ratio of the single-particle 
separation energy to the total binding energy decreases to zero as the number of particles 
increases. The numbers E^"^ /E2 = 16.5 and E^^ /E^'^ = 11.9 obtained from the exact 3- 
body and 4-body results in Eqs. (|3| and (|4) appear to be converging toward the universal 
prediction for large A^ in Eq. ([s]). In Ref. p2j, E^ was explicitly calculated up to A^ = 10 
in lattice effective field theory and found to be consistent with Eq. (|5]). In any real phys- 
ical system, however, the relation ^ can only be valid up to some maximum value of N 
determined by the range of the underlying interaction. When the states become compact 
enough that short-distance properties are probed, the binding energy will no longer be uni- 
versal. In particular, for Lennard- Jones potentials and realistic Helium-Helium potentials, 
effective range effects can be quite large for three and more particles and the universal limit 
is approached only slowly [13j . 

In experiments with cold atoms in a trap, the quasi-2(i limit can be reached by special 
trap geometries. The influence of a trapping potential on ultracold gases in this limit was 
extensively studied by Petrov and collaborators [Ti] - [TB] . Although these works are mainly 
concerned with many-body effects in two-dimensional systems, they have applications for 
few-body aspects as well. 

In this paper, we concentrate on strictly two-dimensional systems of bosons neglecting 
any trapping effects. We calculate three-body observables close to unitarity in the framework 
of an effective field theory for large scattering length. We are especially interested in the 
leading non-universal corrections due to effective range effects. They enter at next-to-leading 
order in the effective field theory. Such effects must be under control for the experimental 
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observation of universal phenomena in 2d. In particular, we calculate the leading non- 
universal corrections to the three-body binding energies, the boson-dimer scattering phase 
shift and effective range parameters, and the three-body recombination rate for finite energy. 



II. METHOD 



In this section, we briefly review the derivation of the three-body equations for d = 2 in 
effective field theory. (See, e.g., Refs. [HI HZ] for more details.) For convenience, we set H = 1 
in the following equations. We include a boson field \l/ and an auxiliary dimer field d in the 
Lagrangian. Since we include effective range effects, the dimer field is dynamical: 

where the dots indicate higher order terms, m is the mass of the particles, rj = ±1, and A 
and g denote the bare coupling constants. The sign rj can be used to tune the sign of the 
effective range term. Negative rj leads to positive values of the effective range r^. In this 
case, the dimer kinetic term has a negative sign and the dimer field is a ghost. We will come 
back to this issue below. Note that three-body interactions enter only at higher orders and 
are not considered in this work. 

The 2d effective range expansion, Eq. ([T]), can also be written in terms of the binding 
wave number k, = ^mE2- 

cot5(A;) = -lnf-V^^(«^' + ^') + (^(^')- (7) 

TT ) 271 

We can deduce the dependence of the binding wave number on the scattering length and 
the effective range from Eqs. ([T]) and ([T]), 



K = --J2W [-2e-^^E-] ^ A/l + 2e-27i5-, (8) 



/ a V 



where the signs are chosen such that k, > 0. The function W is the product logarithm or 
Lambert VT-function. It is defined as the solution to z = we'^, namely W{z) = w. In the 
limit — )■ the expression for E2 reduces to Eq. ^ and k ^ 1.1229/a. 

The Lagrangian in Eq. ^ implies the following Feynman rules: The propagator for a 
boson with energy kg and wave number k is given by i/{kQ — k"^/ (2m) — ie), where k = |k|. 
The bare dimer propagator is i/{ri{kQ — k^ j (4m)) + A) and the boson-dimer vertex coupling 
is given by —ig/2. Because the scattering length is large, boson loops are not suppressed 
and the bare propagator has to be dressed by boson bubbles to all orders. The full dimer 
propagator can be obtained by solving the integral equation in Fig. [2j This leads to the 
expression 

327r f, fp^M — mpn — ^^1 f'^ / o 9 , \ 1 ^ , s 

iD{p^,p) = -i—{\Ti ^ +-{K^ + mpo-pyA)\ , (9) 

where we have already matched g and A to the effective range expansion, Eq. ([T]). The 
wave function renormalization constant is given by the residue of the bound state pole in 
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FIG. 2: Integral equation for the full dimer propagator (thick solid line). The bare dimer propagator 
and the boson propagator are indicated by double and single lines, respectively. 
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FIG. 3: Integral equation for the boson-dimer scattering amplitude. The boson (full dimer) prop- 
agators are indicated by the solid (thick solid) lines. The external lines are amputated. 



the propagator ([9]): 



327r 
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(10) 



The boson-dimer scattering amplitude is given by the integral equation in Fig. [3} which 
iterates the one-boson exchange to all orders. Using the Feynman rules from above and 
projecting onto S'-waves, we obtain [81 ITT]: 



T{p,k;E) 
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m 2 - K2r2 ^(p2 + p _ 
dqqT{q, k; E) 
a/ (p2 + g2 _ mEY — p^q^ 
|g2 _ fYiE — ie 
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K + mE q' 
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(11) 



where k (p) are the relative wave numbers of the incoming (outgoing) boson and dimer in the 
center-of-mass system and E is the total energy. The amplitude T{p, k; E) has simple poles 
at negative energies corresponding to three-body bound states. A more general discussion 
of the analytic properties of few-body scattering amplitudes in 2d is given in Ref . [18] . 
The three-body binding energies are most easily obtained from solving the homogeneous 



version of Eq. (11) for negative energies E 



B{p-E^ 
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a/ (p2 + g2 _^ mE._ 
g2 + mE^ 
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(12) 



In Eqs. (11) and (12) the effective range r"^ is included nonperturbatively in the denom 
inator of the dimer propagator (|9 



Both equations therefore contain some higher-order 
effective range effects but still correspond to next-to-leading order in the effective field the- 
ory expansion. At the next higher order, where terms proportional to (r2)2 enter, there are 
also contributions from the k^ term in the effective range expansion (|l|), (0 which are not 
included here.^ 



^ Note that in three dimensions the /c^ term enters one order higher and the corresponding equation would 
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The integral equations (11) and (12) can be solved in a straightforward way for negative 
effective range [f] = 1). For positive effective range [r] = —1), an unphysical deep bound state 
pole appears in the dimer propagator doj). As the effective range is increased, this pole moves 



to lower energies. Its appearance is related to a violation of the Wigner causality bound 
which constrains the value of the effective range for short-ranged, energy-independent 
interactions. For a detailed discussion of this bound in general dimension d, see Refs. [T91I20] . 
This deep pole appears when we circumvent the Wigner bound by introducing a ghost dimer 
field (?7 = —1). It limits the energy range where our approach is applicable. Identifying the 
position space cutoff in [19] with 1/A, the Wigner bound translates to 



A2 



{[ln(Aa) + l/2]W4} 



(13) 



where A is an ultraviolet cutoff on the integration wave numbers in Eqs. (11) and (12). In 
the limit A — oo, the constraint becomes < . There are at least two strategies to deal 
with this problem: 

1. Expand the full dimer propagator (|9| to linear order in and treat the range per- 



turbatively. 
order. 



This removes the deep pole and includes all terms to next-to-leading 



2. Keep an explicit wave number cutoff A in equations (11) and (12) such that Eq. (13) 



is satisfied. The unphysical pole then has no effect on low-energy observables. 



Both strategies are applicable for wave numbers \k r \ <^ 1. In the following, we make 



use of Eqs. (11) and (12) and use strategy 2 to calculate three-body observables. A brief 



description of the perturbative treatment is given in Appendix Kl 



III. THREE-BODY OBSERVABLES 

In this section, we present our results for the leading non-universal corrections to the 
three-boson binding energies, the boson-dimer scattering phase shifts and effective range 
parameters, and the three-boson recombination rate for finite energy. Since we are mainly 
interested in applications to cold atoms, we will refer to the bosons as atoms in the remainder 
of the paper. 



A. Three-body binding energies 

We start with the effective range corrections to the three-body binding energies in Eq. (|3 



For r^K^ < 0, the energies can straightforwardly be obtained by solving Eq. (12). This is 
similar to the three-dimensional case investigated in [21]. For r^K^ > 0, we have to keep 



track of the Wigner bound. We use an explicit wave number cutoff A and vary A from 1/5 



to 4/5 of the maximum value determined by the position of the unphysical pole in Eq. (12) 



This value agrees within a factor of two with the maximum value given by Eq. (13). The 



be valid to next-to-next-to-leading order. The difference is due to the form of the effective range expansion 
in two dimensions. 
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FIG. 4: (Color online) Three-body binding energies and £'3*^^ in units of h?K^/m vs. the 
two-body effective range r^AC^. The shaded bands are derived with the help of cutoff variation as 
described in the text and provide an error estimate. The crosses are Monte Carlo results for the 
modified KORONA potential from p^[23]. 

dependence of the three-body energies on the cutoff is monotonic. The smallest cutoff results 
in the smallest energy value whereas the largest cutoff results in the largest energy. This 
cutoff variation allows us to check whether the calculation is converged with respect to the 
cutoff and gives an error estimate for our results. We note that one still has to be careful 
about possible artefacts from the iteration of range terms. In the 3d case, it was shown that 
the ultraviolet behavior of the integral-equation kernel is already modified for momenta well 
below 1/r [22]. 

Our results for the three-body binding energies E« and eP function of the effective 
range, r^K^, are summarized in Fig. |4| For negative effective range r^/t^, both three-body 
states become less bound as |r^fi;^| is increased. This behavior is quantitatively similar 
to the dimer state, cf. Eq. ([s]). The binding energies are very sensitive to the effective 
range. For r^K^ = —0.01, we find the values E^^"^ = 1.145(1) for the excited state and 
E^^ = 10.578(1) i?2 for the ground state. For this rather small effective range, the ground 
state energy has already shifted by about 30% while the excited state energy is shifted by 
about 10%. This sensitivity is partially related to the special nature of the effective range 
term in 2d which has units of [length]^. Taking the square root, the leading range correction 
of 10 — 30% for \rK\ = 0.1 looks more natural. Our calculation including the leading non- 
universal corrections suggests that the three-body states eventually cross the atom-dimer 
threshold as the effective range is made more negative. For the excited state this happens 
around r^K^ ^ —0.4, but higher order corrections are expected to be important. If this 
behavior holds true and the effective range could be varied in experiment, the three-body 
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states in 2d might be observable through zero energy scattering resonances similar to Efimov 
states in 3(i. For positive values of the effective range, the central value of our error band 
also corresponds to a less strongly bound system but we can not make a definite prediction. 
Once the effective range effects become appreciable, the errors in our calculation become too 
large. However, the Monte Carlo results for the modified KORONA potential from [131 [23] 
(see the crosses in Fig. |4]) are in good agreement with our band. They also give a diminishing 
energy for larger positive effective range. Note that we only show the two points from [13] 
closest to the unitary limit. All other data points are outside the range of r^K^ displayed in 
Fig-i 

For effective ranges close to zero, the binding energies depend linearly on r^K^. We 
can determine the coefficient of the leading term numerically to about 15% accuracy. For 
r^K^ < 0, we find: 

Ef^/E2 = 16.522688(1) + 28000(5000) r^/t^ + C'(r^K^) , 

eJ^V^2 = 1.2704091(1) + 540(80) r2/t2 + C'(rV). (14) 

For r^/t^ > 0, the coefficient of r^/t^ can not be extracted from our calculation. While the 
values of E^^ and E^^"^ are very insensitive to the cutoff variation at small positive r^K^, 
even the sign of the slope is not well determined. We note that r^K^ can be quite small in 
2d systems even if the effective range is substantially larger that the range of the interaction 
(see, e.g., the explicit example of a circular well potential given in Sec. 6.1 of Ref. [20]). The 



large coefficients in Eq. (14), however, cause significant effective range corrections already 



for Ir^K^I of order 10 ^ to 10 



B. Atom-dimer scattering 

Next, we consider the effective range corrections to elastic atom-dimer scattering. We 
find scattering observables in general to be less sensitive to the unphysical deep poles. To 



obtain the scattering amplitude, we solve Eq. (11) for E = -^k — E2 below the dimer 



breakup threshold with the incoming particles on-shell. The elastic scattering phase shift 
^ad(^) can then be obtained from the scattering amplitude for p = k using 

r(M;i*^-^.)^^/.^^^^^^^. (15) 

In Fig. [sj we show cot 6 ad ^ot different values of r^K^ as a function of the wave number k up 
to the dimer breakup threshold k ^ 1.15^;. For small values of k, cot 6 ad is almost linear in 
Ink but at about half the breakup wave number the behavior becomes more complicated. 
For r^K^ < 0, the non-universal corrections increase cot 6ad compared to the universal result 
while for r^/t^ > the behaviour depends on the value of r^/t^. 

The atom-dimer effective range parameters can be extracted from our results by fitting 
the effective range expansion, Eq. ([T]), to cot 6ad- We have performed fits with different 
orders in the expansion and different truncations of the data sets to estimate the error in 
this extraction. Our results for kqad and (^r^j))^ in dependence of r^/t^ are summarized 
in Fig. [6] The effective range r\j^ comes out positive for all values of r'^K^ considered. For 
|r^K^| < 10~^ the curves are nearly symmetric around r^/t^ = and can be approximated by 

KttAD = 2.614(1) - 4100(500) r\'^ + 0{r^K^) . (16) 
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FIG. 5: (Color online) The elastic atom-dimer scattering phase shift cot Sad for r^K^ = 
0, ±10"^, ±10"^ as a function of the wave number k/n. 
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FIG. 6: (Color online) The atom-dimer scattering length kuad and effective range (urAD)'^ as a 
function of the two-body effective range r^K^. 



Similar to the bound state case, we find large coefficients in the perturbative expansion 
in r^K^. For larger values of r'^K^, the curves are not symmetric anymore. In the unitary 
limit, we find the effective range parameters kqad = 2.614(1) and {KrADY = 4.0(2) . Con- 
verting to units of the scattering length a, our results correspond to a ad = 2.328(1) a and 
In^aAD/a) = 0.845(1). These numbers agree well with the value \n{aAD/(i) = 0.8451 ob- 
tained by Kartavtsev et al. [TT] and are in qualitative agreement with the value qad = 2.95a 
found by Nielsen et al. [6]. For r^K^ of order 0.01, there are again substantial effective 
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range effects. In particular, we find the values naAo = 3.540(1) and (^r^D)^ = T.4(2) for 
r'^K^ = -0.01 and kgad = 2.713(1) and {urAo)^ = 2.9(2) for r^n^ = 0.01. 

We have also calculated the atom-dimer scattering phase shifts and effective range pa- 
rameters using the fully perturbative treatment discussed in Appendix |Aj For sufficiently 
small effective range, the two methods agree. 



C. Three-body recombination 

Finally, we consider three-body recombination into the shallow dimer described by 
Eq. Cold atoms typically also have a large number of deep dimer states. However, 
three-body recombination into the deep dimers is suppressed since the atoms have to ap- 
proach to distances comparable to the size of the deep dimers. This process enters at the 
same order as short-range three-body interactions. It could be calculated by introducing a 
complex three-body parameter [17] . 

The three-body recombination into the shallow dimer is strongly influenced by the be- 
havior of the full dimer propagator, Eq. (|9]). There are three limits in which the propagator 
vanishes: (z) a — )■ oo, (u) E 0, and {iii) E — )■ oo. In these three limits, the three-body 
recombination rate also vanishes. For large but finite scattering length and finite energy, 
however, three-body recombination can take place. 

The rate can be conveniently calculated using the inelastic atom-dimer scattering cross 
section [21] . The integration measure of the three-body phase space in two dimensions using 
hyperspherical variables is given by 

(fpi(fp2(fp3 = ■m^Esm{2a3)dEda3d(pi2dip3^i2(fptot ■ (17) 

From this, the relation between the hyperangular average of the recombination rate at finite 
energy and the inelastic cross section in two dimensions is obtained as 



KiE) = ^ka^Z'\E), (18) 



where k = ^^{E + E2). The inelastic cross section in Eq. (\18n can be obtained by sub- 
tracting the elastic cross section 

^iiiE) = l\h\', (19) 

from the total cross section cr^2'^\E). The latter can be obtained from the optical theorem 
which in 2d is given by ^25j 

^iD^(^) = ^W.(0). (20) 



Using the expression for the elastic scattering amplitude, Eq. (15), the three-body recombi 



nation rate K{E) can be calculated from Eqs. (18), (19), and (20). Our results for K{E) 
as a function of the wave number k are shown in Fig.^ For positive r^/t^, we again make 
use of a cutoff variation in the range of 1/5 to 4/5 of the maximum allowed value. The 
grey band gives our result for r^/t^ = 10~^. For larger values of r^/?^, the width of the band 
increases. As expected, the rate vanishes for large wave numbers and at threshold. Around 
k ^ 3k — 4k, there is a maximum in the recombination rate. The position of the maximum 
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FIG. 7: (Color online) The energy-dependent three-body recombination rate K{E) in units of 
h/{mK'^) as a function of the wave number k/n. The shaded band corresponds to r^K^ = 10~^ and 
is derived with the help of cutoff variation as described in the text. 



is weakly dependent on the value of the effective range. It is governed by the behavior of the 
full dimer propagator as a function of the energy, Eq. ([9]). The maximum is not related in 
a simple way to the energy of the universal three-body states, Eq. ([s]), but our calculation 
establishes an implicit relation between the two. 

In experiments with cold atoms, one typically uses ensembles of atoms in thermal equi- 
librium. The energy-dependent recombination rate K{E) can be converted into an energy 
averaged rate by performing a Boltzmann average as described in Ref. [21]. Taking into 



account the energy dependence of the three-body phase space, Eq. (17), we have 

/;°rfEEe-^/(^^^)ir(E) 



a{T) 



3! dEEe-^I^^BT) 



Our results for a{T) for different values of the effective range are shown in Fig. [8] We do not 
give results for r^K^ > but the shaded band for r^K^ = 10~^ in Fig.[7]translates into an error 
band around r'^n? = here as well. The temperature dependent rate also has a maximum 
at temperatures of the order of 5... 7 times the dimer binding energy. The recombination 
rate at the maximum is very sensitive to the value of the effective range. If the effective 
range is changed from zero to r^K^ = —0.01, the rate at the maximum changes by a factor of 
two. For all values of the effective range considered, however, the recombination rate at the 
maximum remains of order one in natural units h/{mn'^). This suggests that 2d Bose gases 
are stable enough to observe universal few-body phenomena experimentally. The lifetime 
of a 2d Bose gas with large scattering length was previously estimated by Pricoupenko and 
Olshanii [26|. The order of magnitude of our recombination rates is consistent with their 
results. 
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FIG. 8: (Color online) The three-body recombination rate a in units of h/{mHp') in dependence of 
the temperature T in units of h?K^ /{mkB)- 



IV. SUMMARY & CONCLUSIONS 

In this work, we have investigated the three-body properties of identical bosons close 
to the unitary limit in two spatial dimensions. Within an effective field theory for reso- 
nant interactions, we have calculated the leading non-universal corrections which are due to 
the two-body effective range. In particular, we have calculated the leading corrections to 
the three-body binding energies, the atom-dimer scattering phase shift and effective range 
parameters, and the three-body recombination rate at finite energy. 

We have compared our results to previous calculations in the unitary limit where available 
and generally found good agreement. Our calculations show a large sensitivity of three-body 
observables to the effective range. Significant effective range effects can be observed already 
for |r^K^| > 10~^ — 10~^. These corrections are due to large coefficients in the perturbative 



expansion of observables in r^K^ (cf. Eqs. (14), (16)). It would be interesting to understand 



the physics behind these large coefficients. These coefficients could be reduced by an order 
of magnitude by renormalizing to the three-body ground state energy instead of ^2 |27], 
but they would still remain unnaturally large. Our results suggest that the approach to the 
unitary limit in 2d three-body observables is rather slow and effective range corrections play 
an important role even close to the unitary limit. This is in agreement with the results of 
Blume who investigated the universal properties of iV-body droplets using Lennard- Jones 
potentials and realistic Helium potentials [T3] . 

Our calculation of the three-body energies including the leading non-universal corrections 
suggests that the bound states eventually cross the atom-dimer threshold as the effective 
range is made more negative. If this behavior holds true when higher orders are included, 
it opens the possibility to observe three-body states in 2d through a variation of the 2d 
effective range. The states would then appear as zero energy scattering resonances similar 
to Efimov states in 3d. 

Our results are directly applicable to two-dimensional Bose gases with large scattering 



12 



length and imply that effective range effects must be under control in experiments on uni- 
versal properties of 2d Bose gases. Effective field theory provides a powerful tool to calculate 
these corrections and our study provides the first step towards accurate calculations of these 
effects. 

On the experimental side, there has been some progress in the study of universal proper- 
ties of 2d Bose gases. For example, Chin and coworkers have recently studied scale invariance 
and critical behavior near a Berezinsky-Kosterlitz-Thouless phase transition in 2d and ob- 
served universal behavior of the thermodynamic functions [28]. They have also attempted 
to show how the Efimov resonance in three-body recombination shifts when the system is 
tuned towards a dimensionality of two by increasing one of the trapping frequencies in a 3(i 
experiment p9]. 

Interesting few-body properties of 2d systems include universal A^-body states and a geo- 
metric spectrum of A^-body ground states [H]. Moreover, Nishida and Tan have shown that 
a two-species Fermi gas in which one species is confined in 2d or Id while the other is free 
in the three-dimensional space is stable against the Efimov effect and has universal prop- 
erties [30]. More complicated multispecies Fermi gases with similar properties are possible 
as well. They also showed that a purely S'-wave resonance in "id can induce higher partial 
wave resonances in mixed dimensions and pointed out that some of the resonances observed 
in a recent experiment by the Florence group [31] can be interpreted as a P-wave resonance 
in mixed 2d-3d dimensions [32] . Thus, future experiments considering few-body phenomena 
in lower-dimensional ultracold gases will be very interesting. Our calculation of the leading 
non-universal corrections provides a basis for the interpretation of such experiments. 
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Appendix A: Perturbative treatment 



It is also possible to include the effective range in a fully perturbative way. In order 
to achieve this, we expand the boson-dimer scattering amplitude into a leading order piece 



Tq(p, k; E) which satisfies Eq. (11 ) with r = and a correction T2(p, k] E) of order 



,2. 



T(p, k- E) = To{p, k- E) + T2(p, k-E) + .... 



(Al) 



Next 

collect all terms of order r 



we insert this expansion (Al) into Eq. (11), expand the r dependent terms, and 
^ in order to obtain an equation for T2{p, k; E). 



The on-shell 

scattering amplitude at next-to-leading order can then be written as an integral over the 
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leading order amplitude and we finally obtain: 



T{k,k;E) 



1 + 



mr 



2 

2 /"oo 







To{k,k;E) 



< q [ + mE — -g^ 



Toik,q;E) 



In - 777,^; - ie)/K'^) 



(A2) 



where E = 3fc^/(4m) — £'2- The scattering phase shift can be extracted using Eq. (15) as 
before. 
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